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In this paper it is shown how the generating functional for Green's functions 
in relativistic quantum field theory and in thermal field theory can be evaluated 
in terms of a standard quantum mechanical path integral. With this calculational 
approach one avoids the loop-momentum integrals usually encountered in Feynman 
perturbation theory, although with thermal Green's functions, a discrete sum (over 
the winding numbers of paths with respect to the circular imaginary time) must be 
computed. The high-temperature expansion of this sum can be performed for all 
Green's functions at the same time, and is particularly simple for the static case. 
The procedure is illustrated by evaluating the two-point function to one-loop order 
in a (/>g model. 
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I. INTRODUCTION 



Using Schwinger's proper-time representation , the regulated generating functional for 
Green's functions in relativistic quantum field theory and in thermal field theory can be 
formulated in terms of the matrix elements {x\e~^^^\y) and (x|e~^*|?/), respectively, where 
H is generally of the form H = — A)^ + V. Usually these are evaluated using the 
Hamiltonian formalism leading to Feynman integrals for Green's functions. 

On the other hand, in Ref. |^ Schwinger has demonstrated how in particular the one- 
loop effective action can be evaluated by using quantum mechanical methods for the matrix 
elements appearing in the proper-time formulation. In Ref. 0] , it has recently been shown by 
one of the present authors how these matrix elements can be computed by standard quantum 
mechanical path integrals, and that these can be employed to perform perturbation theory 
at arbitrary loop order, using the formalism of operator regularization 0. Independently, 
Strassler motivated by the techniques developed in string theory, has proposed an al- 
ternative formulation based on quantum mechanical path integrals, restricted, however, to 
one-loop order. First-quantized string theory is usually formulated through path integrals 
over fields defined on two-dimensional world-sheets. It therefore seems only natural to expect 
the possibility of rephrasing ordinary field theory in terms of path integrals over functions 
defined on one-dimensional world lines. Indeed, such an approach has already been outlined 
by Polyakov in Ref. |^. The parallel with string theory can be drawn in particular for the 
one-loop effective action in quantum field theory, where there are no internal vertices. Nev- 
ertheless, as shown in Refs. [§,0 in a number of explicit examples, the calculational approach 
based on quantum mechanical path integrals can be exploited also at higher loop order and 
for arbitrary fields. 

In the present paper, we shall show how this technique can be extended to the evaluation 
of Green's functions in the imaginary-time formulation of thermal field theories. Again 
we find that loop-momentum integrals are obviated, however because of the topological 
constraint from the circular imaginary time, an infinite sum corresponding to the winding 
number of paths in imaginary time appears. This sum is equivalent to but not identical 
with the conventional sum over loop Matsubara frequencies. In fact, we shall show that 
the representation as a sum over winding numbers suggests a novel method of evaluating 
the high-temperature expansion of Green's functions. This expansion can be done for all 
Green's functions at the same time, with the coefficients being integrals over the parameters 
corresponding to the insertion of external vertices. This turns out to be comparatively 
simple for the static (but momentum-dependent) case. The procedure is worked out for a 
0^-model, and explicit results are derived for the one-loop two-point function. 
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II. LOOP DIAGRAMS WITHOUT 
LOOP-MOMENTUM INTEGRALS 



If we consider a field theoretical model with an action 



(1) 



(metric signature (— , +,...,+) ) then by introducing a background field [0| /(x) the gen- 
erating functional for one-particle irreducible Green's functions can be formulated in terms 
of the operator 



H 



+ + Xf 



(2) 



where p = —id. 

The one-loop generating functional is given by 



iT^'\f] = lJdx{x\\nH\y), 



(3) 



whereas at higher loop orders {x\H '^\y) is the basic building block. 

In operator regularization H] both In if and are defined through derivatives of H~ 



s-*0 ds 

H-^ = \im^is-l)H- 
s^^ ds 



which is represented as 



H- 



r s) Jo 



dit{it) 



(4a) 
(4b) 

(5) 



In this way the regulated generating functional (which turns out to be finite upon taking 
the limit in s) is expressed entirely in terms of matrix elements {x\ exp —iHt\y); the one-loop 
generating functional in particular is derived from the ^-function M 



T(s) Jo 



dit{it) 



s-l 



according to 



(6) 



(7) 



(Operator regularization at higher loop orders is discussed in Ref. 0.) 

The perturbative evaluation can proceed by using the Schwinger expansion 
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+ . . . 



Jo 



iuHot 



(8) 



and inserting complete sets of states, which involves having to compute loop-momentum 
integrals. 

An alternative for the evaluation of the matrix elements (x| exp —iHt\y) is to use the 
standard quantum mechanical path integral [|l^ 



(a; I exp —i 



-ip-Af + V 



j Dq(r) exp i J dr 

+q{T) ■ A{q{r)) - V{q{r)) 



(9) 



with g(0) = y, q{t) = x. Using a perturbative expansion of the right hand side of Eq. @, 
we can write for the matrix element {x \ exp —iHt\y) with H given by Eq. (|^) 



exp-'- [p^ + Xf]t\ 



Dq{r) E 



N N 



N=0 



i=l 



n / drJiqin)) 



X exp i dr 
Jo 



(10) 



where we have omitted the trivial factor exp —i{m?/2)t. Upon Fourier decomposition of the 
background field, 



(27r)^/2 



d kif{ki)exp-iki ■ qiji), 



(11) 



this matrix element can be expressed in terms of the standard result pO| , pn| for a particle 
subject to an external force 



with 



J Dq{T) expi J dr 

-{x-yf 



^g^(^) - 7iv(r) ■ g(r) 



exp I 



2t 



(27rit)^/2 
1 /•* 

-- / dr-f^^r) ■ [xT + y{t - r)] 
t Jo 

+ 1 f dT'dT"^^{T')-^^{T")G{r\r") 



N 

lN{r) = ^5{t- Ti)ki 

i=l 



2 Jo 



(12a) 



(12b) 
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and G being the one-dimensional Green's function 

1 1 tt' 

G{t,t') = -\T-T'\--iT + T') + — 



(12c) 



By Eq. ( 12a ), the contribution to Eq. (10) that is proportional to when x = y, which 
determines the one-loop A^-point Green's function, is 



-iX 



X 



A^! V2(27r)^/2 

/ rfri . . . rfTAT— — - 
JO [ZTnt 



N ^ N 



(27rzt)^/2 



X exp I 



N 



N 



'^ki ■ X + - ^ hi ■ kjC^Ti, Tj] 

1=1 ^ i,j=l 



(13) 



When substituted into Eq. (|^), the integral over x implements momentum conservation 
ki = 0), and the integration variables Xj, appropriately combined and rescaled, turn out 
to correspond to the usual Feynman parameters. At no stage, however, does one encounter 
any loop-momentum integral. 

This procedure can be applied beyond one-loop order and to models involving spinors 
and vectors An advantage of using this formalism when dealing with a non-Abelian 

vector field theory is that no complicated three or four point vertices need be handled, as 
can be seen from Ref. 1311 . 



III. THERMAL GREEN'S FUNCTIONS 



We now show how this technique can be used to evaluate thermal Green's functions. If 
one were to determine these quantities associated with the model of Eq. ([l|), a starting point 



is the partition function [|T2 



Z = f D(f) exp I / dxQ I ^x 

J periodic L JO J 



X 



period: 
1 

2 



+ (V0)" + m202 



dxo^ 



X 



(14) 



where "periodic" denotes the constraint that (p is periodic in imaginary time xq, 

0(xo,x) = 0(xo + /?,x). (15) 

The regulated generating functional for one-loop Green's functions is again determined 
by a C-function analogous to Eq. (||), 



Vis) Jo 

X J dx dy 5^ (x — y){x\exp —Ht\ 



(16) 
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where the contained in H is contracted with a Euchdean metric. 

The central matrix element to be evaluated and its path integral representation now are 



= (x|exp-- W + \f 



t\ 



DqiT)exp- I 'dr^ [g'(r) + A/(g(r)) 



with g(0) = y, q{t) = x, and with the identification 

qq = qo + n(3, n integer. 



(17) 



(18) 



Paralleling the steps of the previous section, one can again perform a perturbative ex- 
pansion in A and restrict oneself to plane waves 



//5(27r)^/(g(ri)) = exp -i {u^Soin) + ■ q(ri)) , (19) 
where the periodicity condition ( p^Sf ) restricts ci;„. to the discrete set of Matsubara frequencies 

27rn 



(20) 



An A^-point Green's function is thus determined by the expression 



xy 



-X 



N\ \2{27i)(^-^)/^^/p J Jo 



N J 



J Dqir) exp — j dr 



^g^(^) -lN{r)-q{r) 



(21a) 



(21b) 



where 



TV 



i=l 
N 



(21c) 



i=l 



However, because of the topological constraint (|TBp, it is no longer possible to perform a 
shift of variable 



q{r) = q{r)- / rfr'G(r, r')7Tv(r') 



(22) 



with G given by Eq. ( |12c|) to bring the path integral ( pib|) in Gaussian form. Instead of 
a free particle subject to an external force, the quantum mechanical problem described by 
the path integral ( |21b| ) is now that of a particle on a circle (with respect to go)- Its path 
integral treatment has been given in Ref. 11^, and the solution, adapted to our case, reads 
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oo 



p/3 _ V P= 



-)QQ 

z-'S/ " (xo+n/3,x)j/- 

n=— oo 



(23) 



This corresponds to summing over the paths from y to x differing in winding number n 
around the circular imaginary time. 

We consequently need only compute a zero-temperature (/5 = oo) path integral with 
g(0) = y, q(t) = x, and qo{t) = Xq + njS, then perform the sum over n, in order to 
determine the finite-temperature path integral ( |21b|) . With (3 = oo, the shift of variable 
(^) is legitimate. 

Upon substituting Eq. ( [2^ ) into Eq. ( |21b| ) at /? = oo, we find that 



P!^bN] = exp 



^ dr [XT + y(t - r)] . 7iv(r) - ^ drW^r') ■ 77v(r")G(r', r") 

2 JO 



t Jo 



X y Dq{T)exp- dr-q'^ir). 
The remaining path integral is given by 



(24) 



Dq{T)exp- dT-q^{r) = {x\exp-—t\y) 



d^k 



Jk-{x-y)-k^t/2 



(2vr; 



D 



-e 2t 



(25) 



By Eqs. 



and 



1 



we obtain 



-At 



N 



dai . . . daN 



(27rt)^/2 



X exp 



|- [(x - yf + 2n/?(xo - i/o) + n'/^'] /2t 

Af AT AT 

-i{x -y) -Yl ^i^i - W ■'^ki - inPYl ^n,(Ti 

i=l i=l i=l 

1 ^ 1 



(26) 



where we have rescaled = cjit. 

The sum over n can be expressed in terms of the Jacobi function 



so that 



-At 



X exp 



rn + 2i^n 



(27) 



N 



(iai . . .da 



1 



X exp 



(27rt)«/2 

{x-yf ^ ^ ^ 



x^, ^ 



2t 

ip ,xo - yo 



- i^Wiix - y) + y] ■ h + - ki- kjG{ait, ajt) i 
i=i ^ i,j=i 



N 



2ti' t 



i=l 



2TXt 



(28) 
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At one-loop order, only the matrix element with x = y is needed, and integration over 
X again enforces J2i ki = 0. In particular, in D = 6 dimensions, the two-point function with 



ki = k 

xexp 



-k2 and cj. 
1 



ni 



OJr, 



ujn2 is determined by 



dt t 



s-2 



r(s) 647r2/3 Jo JO 
(w2 + k2)(-|ai-(T2| + ((Ti 



daidij2 



0"2 



m 



\ ZTT 



(29) 



The symmetric form of the integrand together with the periodicity of the ^3-function allows 
us to simplify daida2 Jq d{ai — 0-2) in Eq. (pQl). At zero temperature (where 63 ^ 1), 
■u = (Ti — (T2 would exactly correspond to the usual Feynman parameter for combining 
two propagators. The sum defining ^3, however, does not correspond to the conventional 
sum over loop Matsubara frequencies. The latter are recovered after employing Jacobi's 
imaginary transformation |I4| 



9. - 



(30a) 



which in our case reads 



^3 f-^E^*^". 



27Tit 



X 



(30b) 



With Eq. ( P0b| ) , Eq. ( pOf ) coincides with the conventional expression for the finite- 
temperature two-point Green's function where the loop-momentum integration has been 
done, identifying 27m//3 with the usual loop Matsubara frequency. It should be noticed that 
no loop momentum integral was ever encountered in the derivation of Eq. (^); the sum over 
loop frequencies is replaced by the equivalent but conceptually different sum over winding 
numbers in Eq. (p3l).P| 



IV. HIGH-TEMPERATURE EXPANSION 



The standard way to evaluate the finite-temperature Green's functions in the imaginary- 



time formalism |]I2[ is to do the sum over loop Matsubara frequencies first, keeping the loop- 
momentum integrals as such, without combining denominators by Feynman parametrization. 



^The difference of these two sums would have been more conspicuous, had we considered fermionic 
fields: there we would have had an alternating sum over integer winding numbers on the one hand, 
and on the other a sum over Matsubara frequencies determined by half-integers. 
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On the other hand, in our calculational approach, the loop-momentum integrals are circum- 
vented, leaving us with integrals equivalent to Feynman parameter integrals and either with 
the sum over loop Matsubara frequencies or the sum over winding numbers, which are Jacobi 
transforms of each other. 

It has recently been shown in Ref. how to obtain a systematic high-temperature 
expansion of the sum over Matsubara frequencies after the other integrals are all performed. 
In the following we shall demonstrate that the representation as a sum over winding numbers, 
which naturally emerges from the computational scheme considered here, is an interesting 
alternative as concerns the task of a high-temperature expansion. Keeping the integrations 
over the vertex insertion parameters o"j to the end, one can in fact aim at deriving a generic 
form of the high-temperature expansions of all (one-loop) iV-point Green's functions at the 
same time. This will turn out to be feasible and comparatively simple, too, for the static 
case = (kj ^ 0). 

Inserting the matrix element ( PS]) into the (^-function representation of a one-loop A^-point 
Green's function, in the static case we are led to compute the proper-time integral 



2nt 



exp 



--jQNt 



E 



exp 



2t 2^ 



(31) 



where 



Q 



N 



1 ^ 



ki ■ kj (IcTj - cTjl - (cTi 



(To 



2(Ji(jj] 



(32) 



contains the dependence on external momenta and vertex insertion parameters. Integration 
of the latter is postponed for the moment. 

The integral in Eq. ( |3TD is readily performed |]T6[, and the sum naturally splits into a 
zero-temperature contribution (the term with n = 0) and temperature contributions. 



IN 



TOO I T 

J AT -I- J 



N 



{QN/2r~r{s 



V) 



+4E 



N 



n=l 



(33) 



with u = D/2 — N and being a modified Bessel function of integer (half- integer) or- 
der for even (odd) dimension D. Since the temperature contribution is free of ultraviolet 
divergences, we have put the regulating parameter s = there; the otherwise necessary 
differention with respect to s in Eq. (|^) just amounts to dropping the overall factor l/r(s) 
in Ch, Eq. (0). 

The low-temperature expansion of in the form ( p3D can be obtained from the asymp- 
totic formula K^{z) ~ ^7i/2z e~^; the sum over n is, however, slowly converging for high 
temperatures. A series expansion in powers (and logarithm) of f3 can be found by the Mellin 
summation technique [0, leading to 
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c+ioo 



roo I I 

x4 / dn!3-''Q''f rf-''-^K,(l3njQ„) 



(34) 



where c is a real number with c > max(l, 2z/) and ({z) is the Riemann (^-function. Whereas 
the series representation (^3|) is defined only for u > — i, the integral representation (|34D 
allows us to analyticly continue to general u. 

In the form (|34D , can be evaluated by closing the contour to the left and computing 
a sum over residues. There is no contribution from the contour integral over the large arc if 
PVQn ^ 1, which is fulfilled in the limit of high temperature, (3 0. For even dimension 
(integer u) we find 



4 = V(-l)^2^-2^+i5^^-^C(2z/ - 2j)p-^'+'^n^- 



+ 



j=0 

2{-Qn/2) 



r(j + 1 



ln^^-^(^(l) + ^(l + ^)) 



r(2j + i)C(2j + 1) ((3VQ^' 



2i 



(35) 



where if) is the Digamma function. For nonpositive v the first sum does not contribute; for 
z/ < 0, the logarithmic term disappears, too, and the last sum starts at j 
A particular A^-point Green's function is determined by 



I dai . . . da n^n- 
Jo 



For example, the two-point function in six dimensions (z/ = 1) is given by (dropping the 
zero-temperature part) 



nV„ = o,k) = ^£dMa,/,^ 



(36) 



which in the massless case involves only elementary integrations, yielding 



A2 r47r2 



1287r3 



-2 



1 



An 



7 - 



c(2j + i: 



2i 



fy^ ' (2j + l)(2j + 3) V47r 



(37) 



Our results are a direct generalization of the high-temperature expansions obtained in 
Refs. ||r8| , p!9[] for space-time independent quantities such as thermodynamic potentials. In 
fact, the same kind of infinite sums occurs there, the essential difference being that here we 
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still have to evaluate the integrals over vertex insertion (or Feynman) parameters. Indeed, 
the effective potential is contained in our results as the special case kj = 0, where Qn = 
independent of the parameters o"j, rendering the remaining integrals trivial. 

In the nonstatic case, however, where one could derive a generalized version of Eq. (p4|), 
the comparative simplicity of the above derivation is lost. There the condition (3^/QN <^ 1, 
which allowed us to close the integration contour in Eq. (|3^ to obtain Eq. (|35|) is not fulfilled 
because then ^/Q ujn ^ (3. It would be wrong to invoke analytic continuation at this 
point and to substitute iojn — > /cq + ie, with feg a continuous variable independent of /3, 
thereby dropping the contributions from the large arc after all. This supposed transition 
to the real-time formulation is problematic because Eq. has still to be integrated over 



the Feynman parameters implicit m Qn. It has recently been shown by Weldon ||20| that 
although in the imaginary-time formulation Feynman-parametrized Green's functions are 
well-defined, they have to be modified in the real-time theory.0 

Let us emphasize that this presents no problem in principle as our formalism is the 
imaginary-time one, and we have verified equivalence to the usual results above. One has to 
be careful, however, in performing the analytic continuation to real frequencies. The latter 
is straightforward only with the final result at hand, i.e. after all the integrals have been 
evaluated. 

To summarize, we have shown that the recently proposed method of evaluating loop di- 
agrams in field theory by the use of quantum mechanical path integrals can be extended to 
the evaluation of thermal Green's functions in the imaginary-time formalism. This alterna- 
tive approach leads to a representation of the latter in terms of sums over winding numbers 
of paths in place of the usual loop Matsubara frequencies. This different but equivalent rep- 
resentation was shown to also suggest alternative methods in performing high-temperature 
expansion, and we have been able to derive a universal expansion for static one-loop Green's 
functions. 

The extension of our analysis to more interesting field theories is mostly straightforward. 
In Refs. it has been shown how for example Yang-Mills theory can be reformulated 
in terms of quantum mechanical path integrals, if a Feynman-type background field gauge 
is used, leading to a very economical calculational scheme for the evaluation of Green's 
functions. It is not equally straightforward, though, to employ gauges that lead to a more 
complicated kinetic term. Further, with operator regular izat ion it is almost mandatory to 



use a background-covariant gauge fixing p2|. A somewhat more involved case is Green's 



functions with external fermionic lines, because the different boundary conditions of fermions 
and bosons have then to be accommodated in one path integral expression. We intend to 
return to this subject in a forthcoming publication. 



^The uncorrected formulae of Feynman parametrization in the real-time theory as in Refs. |21] 
correspond to an incorrect analytic continuation which has a deceptive analytic behavior at /cq = 
= k, but possesses branch points at unphysical locations, and is singular at infinity [^]. The 
same phenomenon would occur by dropping the contribution from the large arc after closing the 
contour in Eq. (U). 
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